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1 Introduction 



The interplay between quantum and classical aspects in the nuclear collec- 
tive dynamics is a long standing puzzle, similar to the one represented by the 
quantum behaviour of a macroscopic variable [1], and is not yet completely 
understood. 

The phenomenological collective models are based essentially on the quan- 
tization of some simple classical systems (rigid body, liquid drop), with pa- 
rameters obtained by fitting data [2]. Therefore, in this approach the nucleus 
is supposed to be almost frozen, because the model Hamiltonian depends on 
a small number of variables. Also, the Hilbert space constructed by quanti- 
zation is rather artificial, because it does not account for all the observables 
of the many-nucleon system. 

In the microscopical approach, the nuclear dynamics is represented in 
the full Hilbert space by the time- dependent Hartree-Fock (HF) , or Hartree- 
Fock-Bogolyubov (HFB) equations [3]. These equations have static solutions 
corresponding to the ground state, while the excited states are obtained by 
the random phase approximation (RPA), or more generally, using requanti- 
zation methods [4, 5]. A clear sign of collective behaviour, in the sense of the 
phenomenological models, appears when the ground state is non-invariant 
to a continuous symmetry group of the Hamiltonian (spontaneous symme- 
try breaking). This situation may appear for geometrical symmetries, like 
rotation, or dynamical, like particle-number, and the corresponding nuclei 
are known as deformed, respectively superfluid. For these nuclei, the inertial 
parameters of the collective motion may be calculated microscopically using 
the average energy of the ground state, shifted to a "moving frame" by the 
symmetry generators (the "cranking" method). However, the cranking is 
not related strictly to the symmetry breaking and is applied, for instance, to 
obtain the effective mass at fission [6]. 

The choice of the cranking operators is not arbitrary, because they should 
be in one-to-one correspondence with the canonical collective variables, up 
to a unitary change in the Hilbert space representation. This problem of cor- 
respondence is not simple, but a natural solution can be given if the cranking 
is applied at a more fundamental level, to generate both the collective phase 
space and the phenomenological Hamiltonian. The key object necessary for 
this purpose is represented by the shifted ground states, which play only a 
secondary role in the standard cranking calculation. These states are pa- 
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rameterized by the "shift" variables, and may be joined in trial manifolds 
(S), endowed with a classical phase space structure (symplectic form) in- 
duced from the many-body Hilbert space [7] . By construction, S should be 
considered both as collective phase space and as trial manifold for a time- 
dependent variational treatment of the microscopic Hamiltonian. Moreover, 
the Hamilton equations of motion for the collective variables are related kine- 
matically to the evolution of the trial functions produced by the variational 
calculation. Therefore, the "artificial" quantization of the collective motion 
becomes useless, and can be replaced by requantization, to obtain directly 
the microscopic states which correspond to a particular collective motion. 

It is interesting to remark that the standard RPA is included in this for- 
malism as a limit case, reached when S is constructed without any correspon- 
dence to an intuitive model, but taking as generators all possible particle-hole 
(ph) or particle-particle (pp) operators in some finite basis. This nice example 
will be discussed in the next section. Then, Section 3 presents the symplectic 
trial manifolds for deformed and for superfluid systems, and two applications 
to the treatment of the isovector "scissors-like" collective excitations. 



2 Symplectic dynamics on trial manifolds and 



2.1 Symplectic dynamics and requantization 

Let us assume that 7i is the many-body Hilbert space, H is the microscopic 
Hamiltonian, S = {\ip)(X)} is a 2iV-dimensional trial manifold of normalized 
functions, parameterized by the variables X = {x 1 }, i = 1, 2N, and that the 
matrix uj s = [ufJip)], 



is non-singular. Thus, uo defines a symplectic form on S, and the functional 



the RPA 



ufty) = 2klm(d l ^\d^} 




is stationary for the solution X t of the Hamilton equations 



(1) 
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Therefore, the solution of a time-dependent variational calculation within S 
has the form \ip)(X t ), where X t is a trajectory given by Eq. (1). 

It is interesting to note that if the couple H, H corresponds to the quan- 
tum harmonic oscillator, and Sis the manifold of the Glauber coherent states, 
then (1) represents the Hamilton equations for the classical oscillator. 

The procedure of extracting information about the spectrum of H from 
the orbits \ip)(X t ) is called "requantization" , and if the system of Eq. (1) is 
integrable, then the method GIPQ [4] (gauge-invariant periodic quantization) 
can be applied. According to this method, the periodic orbits 7 = {X^}, 

Vl — vl 

should be quantized by a Bohr-Wilson-Sommerfeld (BWS) condition 

[ Tj dt{tl)\id t \^) = 2vrn , n = 0,1,2,.... (2) 
Jo 

This gives the energy spectrum, while the eigenstates corresponding to a 
quantized orbit 7™ should be approximated by the time-average [5] 

|Q 7 ) = J- f Tin dte^mxf) (3) 

with Q t = J* dt? {il>\idt>\i/>) (the "Berry's phase"). 

An extended version of GIPQ includes also the quantization of the invari- 
ant tori [8]. As it will be shown further, the physical applications strongly 
support the requantization by Eq. (3), though its geometrical meaning is not 
yet completely clear. The study of such integral representations ("quantiza- 
tion by membranes") represents an active field of research in mathematical 
physics [9]. 

2.2 The random phase approximation 

We suppose now that the output of a static mean-field calculation (HF or 
HFB) provides the ground state \g) and a set of 2N operators E± a , a — l,N, 
E_ a = E^, so that E_ a \g) = 0. If the single-particle basis contains n states, 
then in the HF case \g) = c\c2cl...c i A \0) is constructed by acting with A 
fermion creation operators c h , h = 1,A, on the particle vacuum |0), and a 
possible choice is E a = c^ p Ch, with 1 < h < A, A < p < n, and N = A{n — A). 
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If there are no "hole" states, and \g) = |0) is the particle (or quasiparticle 
in HFB) vacuum, then we may have E a = c p .c pj , with pi < pj = l,n, and 
N = n{n — l)/2. 

The operators E± a can be used to generate a trial manifold S RPA repre- 
sented by the set of functions 

\1>){Z) = U(Z)\g) , (4) 

where 

U(Z) = eEj^-^E- a ) ( 5 ) 

is an unitary operator and N complex variables. Therefore S is 

parameterized by 2N real variables, denoted x l , % — 1, 2A^, so that for % < N, 
x l are Re(z a ), and for iV < % < 2N, x % are Im(z a ). 

The condition as \g), (or a;* = 0), to be a fixed point for Eq. (1) gives the 
static "mean-field" equations 

{g\[H,E- a ]\g) = , 

automatically fulfilled in the HF or HFB case. 

Let us consider now a small amplitude vibrational periodic orbit 7 around 
\g) with the period T, so that all x k perform harmonic oscillations. This 
means 

with Q = 2n/T, and T, X a , Y" Q unknowns which should be determined from 
the equations of motion. For the orbit 7, 

U{Z2) = exp(e~ int B^ - e iQt B) (6) 

with 

= J2(X a E a -Y:E_ a ) , 

a 

and in the linear approximation Eq. (1) reduces to the RPA-type equation 

(g\[[H,Bi]-mB\E± a ]\g) = . (7) 

The standard particle-hole or particle-particle RPA are recovered when \g) is 
the HF or HFB ground state and E a = c^Ch, respectively E a = c^ p .c^ p .. From 
Eq. (7) one obtains the "normal mode frequency" Q and a one-parameter 
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family of periodic orbits having this frequency. As a parameter we may 
consider for instance the energy 

£={g\U{Z)- 1 HU{Z)\g) . 

However, if S is too large, the amplitudes X, Y increase and Eq. (7) fails in 
approximating Eq. (1). 

Until now the whole discussion was about the classical dynamics on S RPA . 
To establish the connection with the quantum many-body system is necessary 
to requantize the RPA periodic orbits according to Eq. (2). To perform 
this integral we will assume that exists a Hermitian operator W, so that 
[W, B^j = B\ although the explicit form of W is not necessary. If W exists, 
then 

U-Hd t U = n(U~ l WU -W) = Q{e~ int B^ + e int B + [B, 5 f ] + ...) . 

By integrating this sum over a period, the terms linear in B,B^ vanish, and 
the first non- vanishing term is 2n[B, B^]. Thus, the BWS quantization gives 

(g\[B,B%) = n . 

For n = 1 (the first excited state), this coincides with the RPA "normaliza- 
tion" condition, an interesting result proved before using path integrals [10]. 

After quantization, the state associated by Eq. (3) to 7 can be easily ob- 
tained expanding U in powers of B\ B, and retaining only the linear terms. 
The result has the familiar form 

|n> = st|</> , 

but by contrast to the RPA assumption, the excitation operator B' acts on 
the uncorrelated ground state. However, on particular examples it can be 
shown that if the whole expansion of U is considered, then \Q) = B^P RP A\g), 
with Prpa a Hermitian operator which gives an approximate projection of 
\g) on the vacuum \RPA) of B, defined by B\RPA) = 0. Moreover, if \g) 
is symmetry breaking and 7 is a related rotational orbit, then Eq. (3) gives 
P s \g), with P s a symmetry restoring projection operator [5]. 
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3 Isovector excitations in symmetry breaking 
nuclei 



3.1 The two rotor model 

The prediction of the isovector angular rotational oscillations [11] (scissors 
vibrations) in deformed nuclei has been particularly stimulating for the ex- 
perimental research on the nuclear magnetism, leading to the discovery of 
low-lying Ml states. These states have been observed in high resolution 
(e, e') scattering experiments on rare earths [12], /p-shell nuclei [13], and in 
actinides [14]. Their apparent weak excitation in intermediate energy pro- 
ton scattering [15] has supported the orbital character predicted by the two 
rotor model (TRM), but the highly fragmented structure has generated a 
long standing debate about their real origin. On one side were the phe- 
nomenological models supporting the TRM picture, like CSM [16] , or IBA-II 
[17], while on the other were the microscopic RPA or QRPA calculations, 
indicating that the observed Ml excitations are produced by only few quasi- 
particle pairs. Not less important for this debate was the difficulty to decide 
if the states obtained by microscopic calculations correspond or not to angu- 
lar vibrations. Therefore, the problem of finding the appropriate microscopic 
correspondent for a specific collective motion appears to be important. 

This problem can be solved by an RPA calculation based on special trial 
manifolds S rot , generated by cranking, instead of S RPA defined in Section 
2.2. Let us denote by Q x the group of rotations around the X axis, L x the 
orbital angular momentum operator, and by \g) the axially-deformed ground 
state of the microscopic Hamiltonian H. Then, the intrinsic ground state of 
a system rotating around the X-axis with angular momentum C is given by 
the solution \ZJ) of the variational equation 

5(Z\H - uL x \Z) = . 

The set of functions \ZJ) represents a curve in H containining \g) and pa- 
rameterized by the Lagrange multiplier u, or, implicitly, by the angular mo- 
mentum 

£ — \ZJ\L X \ZJ) . 
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The action of Q x moves this curve over a surface in H which contains the 

st ctt GS 

= e-* 1 *^) (8) 

and defines the trial manifold S rot . 

In arbitrary variables {q,p}, the symplectic structure of S rot is given by 
the 2-form oo rot , 

u£ = 2hlm(d q t/j\d p ^) . 

An obvious choice of these variables is q = <p, and p a function of u. When 
this function is the angular momentum, (p = C), then 

= 2nim(d^{(t>,u)\dc^uj)) = d c (Zj\L x \Z u ) = 1 , 

proving that and C are canonical, and S rot is the phase space of the classical 
(plane) rotor. 

In the case of a deformed nucleus, S rot can be constructed separately for 
protons and neutrons, and the trial wave function corresponding to the total 
phase space S™* = S r p ot x S r n ot is 

|V)(^,0„,u; p ,u; n ) =e-<+> I * + +" I M h \Z*) Up \Z n ) Un . (9) 
Let us consider now a schematic nuclear Hamiltonian 

H = E(^o Wj^ - ^-(QisQl + bQi V Ql v ) (10) 

consisting of a spherical oscillator term (ho is the one-body spherical oscillator 
Hamiltonian with frequency c^o = 41 A^ 1 ^ MeV/h), and the quadrupole- 
quadrupole (QQ) interaction, with both isoscalar and isovector components 
(b pa —0.6). Then, for ST°* and H, Eq. (1) takes the form of the Hamilton 
system of equations for two rotors [18], having the cranking moments of 
inertia I p , I n , and interacting by a restoring elastic potential C x (<p p — n ) 2 /2. 
Worth noting, 

C x = 3 Xo (l-b)(Q p ) g (Q n ) g ^ 95 2 A MeV 

appears related to the microscopic QQ interaction (Ql' n = ^/5/167rEp,„(2^ 2 - 
x 2 — y 2 ) Pt n), by contrast to the TRM estimate 2 

C T rm ~ 65 2 A 4/3 MeV , 

2 F. Palumbo, "The scissors mode", in Proc. Int. School Symmetries and semiclassical 
features of nuclear dynamics, Poiana Brasov, Romania, 1986, Springer (1987), p. 230. 
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related to the symmetry energy (5 is the deformation parameter) 3 . Moreover, 
Eq. (2) is identical with the BWS condition for the two-rotor system, giving 
(for n — 1) the quantized angular oscillation amplitudes a T [18], 



1 



*2?il r Ipl n 
- - p, n, I r - 



i \ n ' ±r - i + i ' 

T \ X P ' n 



while the excitation energy is £ x = hQ x = hyC x /I r . 

To apply Eq. (3), \Z) U1 was approximated near \g) using a first order 
perturbative treatment of the cranking term. The "scissors-like" state (not 
normalized) obtained from Eq. (3) has the form 

l^x) = ^(^ + l)[^-«n^]b) (11) 

with D = \6\u . This state gives the B(Ml) strength [18] 

B(Ml) = ^-(g p -g n ) 2 I r D fi 2 N . (12) 

The comparison with the experiment is complicated by the ambiguities of 
separating the orbital and the spin strengths, and defining a reasonable sum 
over fragments. Assuming the dominance of the orbital strength at low en- 
ergy, under 6 MeV in rare earths, and 4 MeV in actinides, the data are well 
reproduced by Eq. (12) with the irrotational moments of inertia. The range 
of the energy S x calculated by taking I r at the rigid or irrotational limits is 
not very large (below 2 MeV) and includes the data throughout all the mass 
regions investigated. 

After normalization, the state \Q X ) becomes the same as the term inde- 
pendent of spin and pairing from the state \ROT), constructed before [19] to 
represent microscopically the scissors modes. However, this is not similar to 
an RPA state, because the "excitation operator" a p L p x — a n L x is Hermitian. 
A quasiboson operator may be obtained if the cranking wave function \ZJ) 
is related to \g) by a unitary transformation. This problem is not easy, but it 
was solved recently [7] in terms of an unitary operator U w which relates the 



3 S = /V45/16^. If (Q p ) g = (Q$) g then 5 = 3 X oCo(Q P ) g /Ku> (c - y/5/4^h/mu> ). In 
C x [18], Xoc 2 Q = 96.3A- 5 / 3 MeV , which yields (Q^ > 3 = (Q P ) g + (Q$) g - 0.18M 5 / 3 fm 2 , 
half the liquid drop estimate AR^S/V^ = 0.36SA 5 / 3 fm 2 , with R = 1.2A 1 / 3 fm. 
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cranked anisotropic oscillator eigenstates to the eigenstates of a s pherical har- 
monic oscillator. Thus, if to is not larger than u> s y/3/2, co s = ^J(^y + uf)/2, 
then we have 

h - -mLU 2 (2z 2 -x 2 - y 2 ) - ul x = U^U' 1 (13) 

with 

U u = e- iXc ° exp(-i Yl d kS2,k) ■ (14) 

k=x,y,z 

In the left-hand side of Eq. (13) the first two terms correspond to an 
anisotropic harmonic oscillator with the frequencies uj x — uj y — u)o\Jl + 25/3, 

oj z = oJo\Jl — 45/3, while 

ho = Y fi^oiblbk + 1/2), = \JmuJo/2h(x k - ip k /mujo) 

k=x,y,z 

h s = Y h^k(blh + 1/2), b\ = ^Jmuj s /2h(x k - ip k /muj s ) 

k=x,y,z 

l X = ih{byb\ ~ b z b\) C X = b\b Z + btby, S 2 , k = l{(b\) 2 ~ (b k f)/A 

and the parameters A, 6 k of \J W are given by 
tan 2A = 2oj/uj s r] sinh# fe = u> s (l — uj k /uj 2 )/2VL k r\ = (cu 2 — uj 2 )/2uj 2 

VL x = uj x Vt 2 yz = (uj s + e y ^ z ) 2 - (uJ s r]/2) 2 

with e y = —e z = u} s i]/2 cos 2A. 

The operators S2,k, k=x,y,z, are the "squeezing" generators, and they 
produce the transition from a spherical to a deformed basis. This transition 
corresponds to the "cranking" of h by the term 6mLo\{2z 2 — x 2 — y 2 )/3, 
proportional to Qq. Therefore, the unitary operator Uq with 9 k as variables 
may be used to generate trial manifolds for the treatment of the quadrupole 
vibrations 4 (S x ~ 2huo). The commutation relations between s 2:k and b\b k 
get closed to the su(l,l) (pssp(l,R)) algebra [20]. 

4 details can be found in Section III of the article M. Grigorescu, N. Carjan, Dissipative 
shape dynamics in the sd shell, Phys. Rev. C 54 706 (1996). 
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The operator c x generates the shift from a static frame to a frame ro- 
tating around the X-axis with the angular velocity u, and it appears as an 
"angle" operator conjugate to l x . Analog operators, c y , c z are associated to l y 
and l z , and by commutation c x , c y , c z and l x , l y , l z generate an su(3) algebra 5 . 
Similarly, S2, x , $2,y, $2,z and l x , l y , l z generate by commutation a gl(3,R) alge- 
bra. This rather complicated set of algebras is included in the symplectic Lie 
algebra sp(3,R) [21]. 

In the many-fermion case the one-body operators h,ho,h s ,lk,Ck,Sk be- 
come particle-hole operators, denoted H, H , H s , Lk,Ck, Sk, and may be used 
to write \Z) U1 as 

|Z) w = e- iAC -exp(-i e ^2,k) \9s) ■ 

k=x,y,z 

This form is especially suited to study large amplitude vibrations, but cor- 
rections to Eq. (11) appear already in the linear approximation [7]. If the 
dependence of 6k on u is neglected, then during the scissors vibration each 
\Z) UJ in Eq. (9) changes in time only due to the factor of C x , and Eq. (3) 
gives \n x ) = St with [7] 



B\ = L[ a P L P a n L n_ im x (aV C v a n Cl)\ . (15) 



Unexpectedly, though a p ' n are the same angular amplitudes as in Eq. (11), 
£?t is normalized in the RPA sense, {g\[B x , B^\\g) = 1. In fact, if the 
particle-hole excitations between different oscillator shells are neglected, then 
sgn(i])ihC x \g) ~ —L x \g), and B x \g) reduces to Eq. (11). 

The operator B x was obtained recently also by the canonical quantization 
of the TRM in relative coordinates [22], and it was proved to support the 
interpretation of all low-lying orbital 1 + excitations as a scissors mode. 

It is interesting to note that if the dependence of 9^ on uo is taken into 
account 6 , then B x will contain operators from gl(3,R). Such terms are related 
to the excitation operator proposed long time ago by Hilton [23] . Instead of 
constructing B % within the "angular momentum & shift" algebra, su(3), he 
proposed a combination within the "angular momentum & squeezing" alge- 
bra, gl(3,R), between l y and the "shear" generator zp x + xp z ~ i(b\b x — b z b x ). 



presented e.g. in math-ph/0007033. 

6 M. Grigorescu, D. Rompf, W. Scheid, Dynamical effects of deformation in the coupled 
two-rotor system, Phys. Rev. C 57 1218 (1998). 
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3.2 The isovector Josephson oscillations 

The combined effect of the proton-neutron interaction and breaking of the 
translational or rotational symmetries is related to the giant dipole reso- 
nance [24] or to the "scissors modes", respectively. However, there is one 
more important symmetry breaking in nuclei for which such type of isovec- 
tor collective motion was not yet observed, and this appears when the nuclei 
are " superfluid" . 

The ground state of a superfluid system accounts for the pp correlations 
produced by the pairing interaction, and is well approximated by a BCS 
function. For a single j-shell, the pairing Hamiltonian and the BCS function 
are 

H = eN-jP^P , \BCS)(<p,p) = e<"«-<' p >\0) , (16) 

where z = pe~ lip , p is the BCS "gauge" angle, iV = J2m=-j c m c m is the 
particle-number operator, and 

p j = t (-ir m 4ci m 

m=-j 

is the pair creation operator. 

The angle p is not a constant, and a superfluid system in the ground 
state performs a free gauge rotation with the angular velocity ip = 2ep/h, 
twice the Fermi frequency. In a nucleus, the proton and neutron systems 
are not isolated, but change particles until the Fermi energies e p F) e F become 
equal. Thus, we may see this as an indication about the existence of a phe- 
nomenological "gauge restoring interaction", which tends to fix the relative 
gauge angle p p — p n at a constant value. If this is true, and there is an 
interaction between pairs of protons and neutrons, then oscillations of the 
protons against neutrons in the BCS gauge space should appear. 

A Josephson-like proton-neutron interaction 7 

H pn = - a -(Pip n +pip p ) 

may be related to the isospin symmetry breaking mean-field for a four particle 
interaction [25]. The problem of the mean- field created by the four particle 

Considered also in M. Gerceklioglu, A model for the doublets of the K 7r =0 + states in 
deformed nuclei, Acta Phys. Slovaca 52 161 (2002). 
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interaction is not new, but previously [26] the main interest was for terms 
~ PpPn, assumed to represent a clusters, while terms like PpP n were ne- 
glected. Because H pn does not commute with the isospin T = (N n — N p )/2, 
it produces an " isorotational" term k a (N — Z) 2 in the total energy. This 
means a term in the symmetry energy kw{N — Z) 2 , (kw = 28/ A MeV), from 
the Weiszacker mass formula with a "dynamical" origin, beside the "kine- 
matic" one determined by the Pauli principle [27]. Thus, a can be obtained 
by a fit of the symmetry energy produced by the Hamiltonian 

H = Hq + Hq + H pn (17) 

in a single j-shell. Considering the case of le/3/2 nuclei, a was estimated 
to be Rj 2.1 1 A MeV [25]. However, an approximation of j-shells with high 
degeneracies suggests a value about ten times larger [28]. 

The interaction H pn contributes to the symmetry energy in all nuclei, but 
in the superfluid nuclei it produces also a restoring potential C a ((p p — (p n ) 2 /2 
for the BCS angles. This potential can be related to H pn by a treatment 
similar to the one applied to the QQ interaction responsible for the orbital 
scissors modes. Indeed, the BCS functions define symplectic manifolds S BCS 
which can be parameterized by f and p, or by the canonical variables ip and 

p= (BCS\N\BCS)/2 = (j + 1/2) sin 2 2p . 

In these variables 

uj^ s = 2hIm(d v BCS(ip, p)\d p BCS(ip, p)) = hd p ((BCS\N/2\BCS)) = h . 

For a proton-neutron system the trial manifold will be represented by the 
product Sp CS x Sn° s ' , and the collective motion determined for H by Eq. 
(1) shows the occurrence of isovector "gauge- angles" vibrations [25]. For a 
half-filled shell, the fixed point in (1) corresponds to the ground state of H, 

|^ )=e f(Pt + iVt-P-iV) |0) (lg) 

and the gap parameter A = G(g\P^\g) /2 is the same both for protons and 
neutrons. The restoring potential has the constant C a = 2a(A/G) 2 , and the 
excitation energy for the isovector oscillations is E a = hQ a = A^/2ky/C a [25]. 
The excitation operator defined by Eq. (3) has the form 

Bl = rTjri^[To - -P-rf + N)] (19) 
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and by accident, is the same as the one provided by the standard QRPA [29]. 
However, the RPA vacuum defined by B a \RPA) = exists only if j + 1/2 is 
an even integer [28] . 

The states generated by are isovector monopoles, and correspond to 
Josephson oscillations between the proton and neutron superfluids 8 . Such os- 
cillations might be excited by the Coulomb interaction in electron scattering 
[28] , or by the current of pairs between the two superfluids produced in pion 
double charge-exchange (DCX) reactions. The importance of Josephson-type 
correlations in DCX reactions was proved first by shell - model calculations 
[30], suggesting that the "scissors modes" in gauge space discussed here are 
worth of experimental investigation. 



4 Summary 

In this lecture I presented a microscopical approach to the collective motion, 
based essentially on the time-dependent variational principle and GIPQ re- 
quantization (Section 2.1), but which is peculiar by the choice of the trial 
functions. The trial manifolds are supposed to have the phase-space (sym- 
plectic) structure of a specific collective model, and for the symmetry break- 
ing nuclei are constructed using the cranking procedure [7]. The Hamilton 
equations of motion (Eq. (1)) appear by a constrained variational calculation 
in the Hilbert space, rather than by semiclassical approximations (h — > 0). 
This formalism includes the standard RPA or QRPA (Section 2.2), and was 
applied with success to the treatment of the collective isovector excitations 
(Section 3). It solves the problems of the inertial parameters, restoring force 
constants, and of the microscopic analog for a particular collective motion 
(here the scissors vibrations). Moreover, isovector vibrations in the BCS 
gauge space of the superfluid nuclei were predicted. The variational formu- 
lation is appropriate to account for the coupling between a quantum system 
and a thermal environment [31]. Therefore, the study of giant resonances 
using a Langevin form of Eq. (1), with noise and memory- friction terms in 
the right-hand side, appears highly interesting. 

8 At <7<Gan odd-even effect in the total number of pairs can appear (M. Grigorescu, 
Low-lying excitations in superconducting bilayer systems, cond-mat/9904242, or High-Tc 
Update 13 No. 10, May 15 (1999), or Can. J. Phys. 78 119 (2000)). 
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